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Abstract: The chiral-odd distribution function hi is discussed in the framework of the 
renormalon approach. Using a bag- model calculation for hi, we calculate its intrinsic 
uncertainty due to renormalon poles. The result is given as a function of Bjorken-x as well 
as for the first moments separately. We also extract the perturbative corrections to the 
first moments of hi. 
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1 Introduction 



While in totally inclusive deep inelastic scattering (DIS) the quark chirality is conserved 
up to terms proportional to the quark masses, this is not the case in the Drell-Yan process. 
Here a quark-antiquark-pair is annihilated to a virtual photon, so that in the cross-section 
the quarks originating and ending in the same nucleon may carry different chirality. This 
gives raise to chiral odd distribution functions, which first appeared in the discussion of 
the transverse polarized Drell-Yan process [0]. This chiral odd distribution function is 
defined by a twist-2 operator and is called the transversity distribution hi. Unlike the 
helicity asymmetry gi, hi has no partonic interpretation in the chiral basis. Changing to 
the transversity basis gi looses its partonic interpretation and hi gets one instead. It 
is interpreted as the probability to find a quark in a transversely polarized nucleon in an 
eigenstate of the transverse Pauli-Lubanski vector with eigenvalue +1/2, minus the same 
with eigenvalue —1/2. 

Some experiments are planned to measure hi in the near future, especially at RHIC 
(BNL) and possibly by the COMPASS experiment at CERN — for a general review of 
the possibilities for measuring the transversity distribution see ^ — , so that there is need 
for theoretical predictions for hi. The anomalous dimension which we will also consider 
in this contribution was calculated in Q]. In the nonrelativistic quark models gi and hi 
are identical. The positivity of parton probabilities implies the inequality |/ii(a;)| < /i(a;) 
and the Soffer inequality |^ 2\hi{x)\ < fi{x) + gi{x). A bag model calculation predict 
|^i(^)| > which may be correct in general 0. For medium large Bjorken-x there 

exists a QCD-sum rule calculation 0, predicting a much smaller value for hi than the 
already mentioned bag model calculation. 

We wish to contribute to this theoretical discussion by calculating an intrinsic uncer- 
tainty of the perturbative series due to IR-renormalon poles. This uncertainty is a system- 
atic error for the experimental measurement of hi. In addition we discuss the anomalous 
dimension in the one loop approximaion and extract the correction to the first moments 
of the transversity distribution. 



where x = — g^/(2p • q) is the Bjorken variable, ra'^ is a light cone vector with n"^ = 
of dimension (mass)'^, p is the proton momentum and s is the proton spin, = m^, 
= — 1 and p ■ s = 0. p ■ n = 1 and the transverse part of the spin vector is defined by 
the decomposition = (s ■ n)p^ + {s ■ p)n^ + s^^. The contraction of this expression with 
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the light cone vector n'^ gives: 

^ J ^e'y+''{p,s\ip{0)a^^n''-f5ip{y+n)\p,s) = g^)s_L^ - g^)m^?2^(s ■ n) . (2) 

It is possible to relate the transversity distribution hi to the imaginary part of the forward 
scattering amplitude on the leading twist-level 0]: 

T, = '-Jd'ye"^^{p,s\T{j,,{y)js{0)+js{y)Mm\p,s) . (3) 

Here j^^iy) = 4'{y)l^il54'{y) is a axial-vector current, js{y) = 4'{y)i^{y) is a scalar cur- 
rent and T denotes the time-ordered product. Summation over flavor indices is assumed. 
Equivalently one may use a time-ordered product of a vector-current jfj.{y) = 4'{y)lfj.i^{y) 
and a pseudoscalar current js^/) = 4'{y)l54'{y)- 

To prove the relation between this definition and the operator-definition Eq. (|l]) of the 
transversity distribution, one has to decompose the forward scattering amplitude into its 
Lorentz-structures using the conservation of either the axial-vector or the vector current. 
The conservation of the axial- vector current is correct only for the flavor nonsinglet current, 
so that the proof remains correct up to an order a^-correction only. As the vector-current 
is conserved independently of the flavor combination under consideration we prefer the 
definition 

T,= '-Jd^ye'''{p,s\Tij,iy)UO)+Uy)j^iO))\p,s) . (4) 

To find the realation to hi the current product has to be expanded collecting the terms with 
one incoming quark, one outgoing quark and one quark propagator S{y) = i ^ A{y), where 
A is the Pauli- Jordan function. The Pauli- Jordan function is expanded on the light cone 
and only the leading term is taken into account. The imaginary part of the s-channel-term 
of the forward scattering amplitude takes the form: 

ImT^ = dSe''y{p,s\ : ^{y)a,,y''t^,m ■ \p,s)S'{y') 

+ ^ I d%e''y{p, s\ : m^^uV'^l^Hy) ■■ b, s)5'iy') . (5) 
Ztt J 

Here 0"^,^ = | [7/i,7z/] and S'{y'^) = {d/d{y'^)) 6{y'^). Substituting light cone variables y- = 
yo — y?j and 2?/+ = yo + 1/3, integrating by parts, using translation invariance, and putting 
this expression on the light cone by choosing y such that = 0, one gets 

ImT^ = -l|^e*^+^(p,s|:^(0V^,n^75^(y+n):|p,.) . (6) 

This result has the form of Eq. (|^), so that one gets a relation between the s-channel-part of 
the forward scattering amplitude and the transversity distribution by collecting the terms 
proportional to s^^\ 

hi{x)s^^ = - ImT^I . (7) 
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2.2 Light Cone Expansion and Moments 

The light cone expansion of the forward scattering amphtude may be writen as: 

T^^l,^ = 2s±^ Y.Cm[^,as] ^-+1 + highertwist (8) 

where u = 1/x. are the reduced matrix elements of the local twist-2 operators relevant 
for hi: 

{p, s\^a^^H-ir^iD^' ■ ■ ■ iD^'^^ip\p, s) = 2Am{s^p^^ - p^s^'')p^^ ■ ■■p'^^^ij - traces . (9) 

The brackets denote total symmetrization of all included indices. The symmetrization 
and the subtraction of the traces are necessary to extract the leading twist part of the 
matrix element. The moments of hi have to be defined as = Cm Am to get the general 
expression: 

2sx^M„ = - /'rfxx"(JmT^(x) + (-r/mT^(-x))l . (10) 

7T Jo 

Using Eq. (0) the moments become: 

Mn = l f' dxX^ihi{x) + i-rhi{~x)) (11) 

/ Jo 

for n = 0, 1, . . . This result coincides with the expression found by Jaffe and Ji 0. 

The operators in Eq. have a well defined behaviour under charge conjugation: they 
are C-odd for even n and C-even for odd n. This change in sign corresponds to the relative 
sign of hi{x) and hi{—x) in the moments Eq. (|TTp. To obtain the correct relation between 
the antiquark and the quark transversity distribution one may look at the first moment: 

Mo = I f'dx (hiix) + hi{-x)) . (12) 
/ Jo 

As this expression is odd under charge conjugation, the antiquark transversity distribution 
should be hi{x) = —hi{—x). It follows that the contributions of sea quarks cancel for 
even moments in general and only the valence quarks contribute, while the sea quark 
contributions add for odd moments. 

3 Renomalon ambiguity of hi 

We are calculating the IR-renormalon contribution P] to the structure function defined 
by Eq. (^) and (0). As its twist-2 part is related to the transversity distribution hi, the 
IR-renormalon contribution may be interpreted as intrinsic ambiguity of the perturbative 
expansion to hi. Perturbative series in QCD are asymptotic ones, implying that the 
radiative corrections of higher and higher orders get smaller only up to a finite order 
mo and diverge for larger orders in a^. The uncertainty of the series is of the order of 
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the smallest contribution and has the form of a power correction ||^, so that the QCD 
perturbative series on the lowest twist level may be written as (a^ = 05/(47?) ): 

a(Q^ = /.^a,) = + + + • (13) 

Note that the 1/Q term is not present in the above expression. It was shown for Drell- 
Yan on the one gluon exchange level that these corrections are cancelled by higher order 
perturbative contributions JlOl and that the 1/Q^-term is the leading power correction. The 



uncertainty may be determined by the exact calculation of the perturbative corrections up 
to the order mo(Q^), which is a very demanding procedure. 

Instead we calculate the forward scattering amplitude Eq. (^ on the one gluon exchange 
level, using a Borel transformed effective gluon propagator O]: 



BvaAa.DfM]{u) = . (14) 



C corrects for the renormalization scheme dependence (C = — | for MS-scheme), fi is the 
renormalization scale, and u is the Borel parameter. The effective gluon propagator is 
constructed by replacing the coupling by the running coupling constant, that is by a 
resummation of all quark- and gluon-loop insertions in one gluon-propagator. In the first 
order of this propagator leads to exact results. Looking at higher order corrections the 
restriction on one exchanged effective gluon corresponds to the large iVj-limit [0, where 
Nj denotes the number of quark flavors. The next-to-leading A'^y-terms are approximated 
by naive-nonabelianization (NNA) ||13[. This corresponds to the replacement of the one 
loop QED-beta- function by the QCD-beta-function /3o = 11 — |A^/ or equivalently to 
Nf — > Nf — y. The quality of this approximation has been checked |]T^ |T5| for the 
unpolarized structure functions F2 and Fl and for the polarized structure function gi by 
comparing the NNA perturbative coefficients with the known exact ones. This comparison 
gave very reasonable results for Fi and (71, while the results for F2 are less convincing. 

Formally, asymptotic freedom is destroyed in the large A^j-limit. One should recognize 
that the large A'j-limit is used to select graphs and has to be understood as a definition 
of an approximation procedure. At the end Nf will be set to 4, so that /5o stays in an 
asymptotic free region. This procedure is technically analogous to the use of the large 



A'c-limit |jT6[, even if the physical content is different. 

In the Borel plane the ambiguity of the truncated perturbative series in Eq. (|IB|) is re- 
flected in IR-renormalon poles, which hinder an unambigous inverse Borel transformation. 
This ambiguity of the perturbative series may be interpreted as twist-4 contribution of 
the structure function defined by Eq. (|^) and (^. Such estimations of higher twist correc- 
tions gave very reasonable results [1T7| , |l^ , 14, |l^. However the transversity distribution 
is defined by pure twist-2 operators and the relation of Eq. (^ to the operator definition 
of the transversity distribution Eq. (|l|) was shown for the twist-2 part of the structure 
function only. This means that the IR-renormalon ambiguity of the structure function has 
to be interpreted as intrinsic uncertainty of the perturbative expansion of the transversity 
distribution hi. 
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So let us calculate the Borel-transformed s-channel forward scattering amplitude in 
Eq. (^) on the one gluon exchange level using the effective gluon propagator Eq. (|14D. The 
result is a series in u = 1/x which has to be compared with Eq. 



oo 



^0 V 




m+1 



+ higher twist (15) 



where the reduced matrix element was determined at the tree level to be y4„ = 1 and a 
factor 2 is missing because the exchange graphs are not included in this expression due to 
the restriction on the s-channel contribution. We obtain for the Borel transformed Wilson 
coefficient: 




Til + s + n) 



E 



T{s + k) 



+ 



+ 



1 - s 
1 



i + s r(i + s)r(i + n) ' ^tl r(i + s)r(i + k)i + s + k_ 

^ _ 2V{1 + s + n) ^ " T{s + k) 4(1 + fc) 

^ ^1 r(l + s)r(l + A;)l + s + fc 



i + s r(i + s)r(i + n; 

1 r(l-|-s + n) 

m ^ r(i + s)r(i + n) ' r(i + s)r(i + k)i + s + k_ 



E 



T{s + k) 



(16) 



0,1,2, 



where s = PqU replaces the Borel parameter u. We find IR-renormalon poles at s 
as usual for DIS. This is not a general statement as in the case of e"'"e~-fragmentation 
one obtains an infinite sum of poles with even powers of 1/Q [ll^. The Wilson coefficient 
has still to be renormalized and the coefficient of the 1 / s-pole should give the anomalous 
dimension as will be discussed in the next section. 

In the Borel-plane the ambiguity C^^ of the perturbative series in Eq. ( [I3l) can be 
rediscovered as ambiguity of the inverse Borel transformation due to the IR-renormalon 
poles or in other words as the imaginary part of the Laplace integral: 



'Ale 



— Im 



ds e-^/(*"^)i3 1 



C-n 



(17) 



We get 



2Cf 
Cf 



3 + n - 



l+n 2+n 




±^ |3 + -(l + n)(4 + n) 



l+n 2 + n 3 + n k=i^ 




The signs of these twist-2 uncertainty terms remain undetermined, because it is not clear 
in which way the pole should be circumvented in the Laplace integral. 

The ambiguity of the Laplace integral is an intrinsic uncertainty of the twist-2 transver- 
sity distribution. The ratio of the moments of this uncertainty and the moments of hi 
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is expanded up to the order ag/Q 



2 



^(1) A^e-g I ( as _l\ 



Mn . (19) 



Here the truncated perturbative series Eq. ( pTSD was used for the moments of hi. The lowest 
order twist-2 perturbative coefficient B^^^ is determined by the tree-graph and is 1. As 
an illustration we will insert for M„ a theoretical model prediction for the transversity 
distribution. 

lYiom. Eq. (0) and ( ]T9| ) we can calculate the IR-renormalon uncertainty f/„ for each 
moment, where 

= 4°)(l±?7„) + 0(a,) , (20) 



where A^^^ denotes the leading contribution to the moments. At = 4 GeV^ and with 
= 200 MeV, A^^ = 4, = | and /5o = 11 - |A^/ we find: 

nl 



dx (hi{x) - hi{x)) = 4°^ + C(a,) , 

£ dxx(hi{x) + hi{x)) = 4°^(1± 0.0057) + C>(a,) , 

J dxx"^ (^hi{x) -hi{x)) = Af\l±0.0U) + O{as) , 

dxx^ (^hi{x) + hi{x)) = 4°^(1 ± 0.024) + 0(a,) , 

dxx^(^hi{x)-hi{x)) = Af\l± 0.036) + (as) . (21) 

We get no IR-renormalon uncertainty for the first moment. The IR-renormalon uncertainty 
becomes larger for higher moments, so that we expect larger ambiguities of the transversity 
distribution in the region of larger Bjorken-x. For higher moments the contribution of hi 
can be considered as marginal, so that the above ambiguities should remain approximately 
correct for dxx"'hi(x) with n > 2. This is of course not the case for the first and the 
second moment {n = 1). A rough estimate gives rise to a sea-quark effect of the same 
order as the calculated uncertainty. 

/^From the moments in Eq. (|l^) the valence quark transversity distribution hi — hi may 
be reconstructed as a function of Bjorken-x. The result is a convolution integral: 

h{^{x,Q')-h{\x,Q') = ^r7^^'H^){^i(^)-^i(^)} (22) 
where C^'^^x) is defined by 

CW = / dxx''C^''\x) (23) 
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Figure 1: The bag model calculation for hi — hi (full line) and the IR-renormalon ambiguity 
evaluated at = 4GeV^ multiplied by a factor of 10 (dotted line). There are two changes in 
sign at a; ~ 0.25 and x ~ 0.4. The ambiguity was added and subtracted from the bag model 
calculation (dashed lines), (Ajg^ = 200 MeV, and Nf = 4). 

and we find for the two first IR-renormalon ambiguities: 

C7«(x) = ±^|— + 35(x-l)-y(x-l)-2a;-2l 
Po i (1 - x)+ J 

c^'\x) = ±^ { (13^ + - 1) - ^'(^ - 1) + - 1) - - X - 1| . 

(24) 

The IR-renormalon ambiguities shown in Fig.[^ are calculated using the bag model 
calculation of for hi{x). The large A^y-limit is most reliable in the region of medium large 
Bjorken-x. For small x the neglection of multiple gluon exchange is no longer justified. 
On the other hand the influence of the hadronic spectrum makes a pure perturbative 
calculation insufficient for large x. In the region of best accuracy (0.2 < a; < 0.45) the 
uncertainty does not become bigger than 1% (see Fig.^). One can expect a sizeable IR- 
renormalon uncertainty of up to 10% for 0.5 < a; < 0.6. 
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Figure 2: The relative magnitude of the IR-renormalon ambiguity with respect to the bag 
model calculation [§ for hi - hi {Aj^ = 200 MeV, = 4GeV^ and Nf = 4). 



a. 



corrections to the moments 



It is possible to extract the anomalous dimension and the corrections to the moments of 
the nonsinglet transversity distribution from Eq. (pIBI). This expression has to be renormal- 
ized so that the Borel transformed Wilson coefficient should be regularized simultaneously 
analytically — like it was done in the previous section — and dimensionally. We get for 
this generalized expression {d = 4 — e): 



B 



Cn 



E 

fc=0 



rri 



+ A;)r (^ 



-4V\^ r(2-|)r(-|- 

^ ) T{3-e-s)T{l + s) 



r(l + A;)r(2 + s + A;) 



s + n + 



T(l+n) 



+ es{2 



s + -j (6{l-s) + {A-e){l-e))+Ak{2 
r(l + s + n+|)r(l + s + n)| 



(25) 



' T{3 + s + n)T{l + n) 

It is important to take care of a possible violation of the Ward identities due to the analytic 
continuation of 75 to ^-dimensions 
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The coefficients of the perturbative series can be reconstructed from the Borel trans- 
formed Wilson coefficient by an appropriate number of derivations with respect to s at the 
origin s = of the Borel plane: 



^(fc+i) 




(26) 



s=0 



This may be verified starting with Eq. ( [T^ ) and using 



(27) 



So the ffist order as correction is determined by setting s = in Eq. (^). We find 
the quark-quark anomalous dimension of the nonsinglet operator as the coefficient of the 
divergence at d = 4: 

This anomalous dimension is identical with the coefficient of the 1/s-term in Eq. ([T6| ) 
in the pure analytical regularization, as it should be. In momentum space we get the 
corresponding splitting function 



Prix) 



Cf_ 
47r2 



[l-x) 



(29) 



which is to be compared with the splitting function for hi as found in the literature 



P 



X] 



Cf 



1 



(1 



X] 




(30) 



Both expressions differ by |(5(x — 1). This term is not really calculated in the conventional 
approach but rather fixed by the requirement that the transversity flip probability for 
quarks emitting a zero momentum gluon has to vanish]^ In addition the (5-function part 
should be the same as for the polarized structure function gi. Following this standard 
argumentation the corresponding 5{x — l)-term has to be added in Eq. (^). However, 
such additional arguments should not be necessary in a direct calculation of the anomalous 
dimension. We got the right 5-function in the case of gi without any additional arguments 



in an analogous calculation ||T5[ . It is clear from a comparison of these two calculations that 
the behaviour of the anomalous dimensions of gi and /ii at x — > 1 may only be indentical, if 
the vertex-correction with scalar- or vector- vertex contributes in the same way for x — > 1. 
This is obviously not the case, so that the difference of Eq. ( 129|) and (|30|) at x = 1 appears 
necessarily. Considering the general arguments leading to Eq. ( ppj ) as correct, we conclude 
that the equivalence of the twist-2 part of the used operator in Eq. (|^) to hi leading to 
Eq. (^91) is problematic at x = 1. One should emphasize, however, that the renormalon 
estimation in the last section is not affected of this question. 

^We thank X. Artru for helpful comments. 
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Unlike in the case of gi the anomalous dimension does not vanish for any n, so 
that at one loop level a renormalization is already necessary for all moments. We choose 
the MS-renormalization scheme. After renormalization of the Ward-identity with 

Z, = i-CF — + 0{a'',) (31) 

TT 

we get for the first moments of the structure function defined by Eq. (^: 

J dx (hi{x) -hi{x)) = l + 0[af) 
[\xx(hi(x) + hi(x)) = l + lcF— + 0(a^,) 

Jo ^ ' 8 TT ^ ^ 

dxx^ (h{x)-hi{x)) = l + ^CF— + 0(a^) 
Jo ^ ' 8 vr ^ ^ 



dxx^(h^[x)^hx[x)] = l + ^CF— + 0(al) 
\ I 48 vr ^ ^ 



/ dxxHh^{x)-h^{x)] = l + —CF— + 0(al) . (32) 
Jo ^ ' 8U vr ^ ^ 

As the used NNA- approximation is exact on the one-loop-level, these first order perturba- 
tive corrections are exact results. 



5 Conclusions 

We calculated the IR-renormalon ambiguity for the valence quark transversity distribution. 
The ambiguity is smaller than 1% in the region of best accuracy of the NNA- approximation. 
For X ~ 0.6 this systematic error becomes important reaching about 10%. Thus the 
renormalon uncertainty should be taken into accout when interpreting measurements of 
h\ in this region of x. In addition we presented a calculation of the anomalous dimension 
using the operator in Eq. (^) and of the corresponding one loop perturbative corrections 
to the first moments. To our knowledge we calculated the perturbative corrections for the 
first time. 

Acknowledgements. We thank L. Szymanowski and S. Hofmann for helpful discus- 
sions. This work has been supported by BMBF. A.S. thanks also the MPI fiir Kernphysik 
in Heidelberg for support. 
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